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I compute the effect on the power spectrum of tracers of the large-scale mass-density field (e.g., 
galaxies) of primordial non-Gaussianity of the form $ = <£ + /nl (4> 2 — (4> 2 )) + <?nl<^ 3 + where $ is 
proportional to the initial potential fluctuations and <f> is a Gaussian field, using beyond-linear-order 
perturbation theory. I find that the need to eliminate large higher-order corrections necessitates the 
addition of a new term to the bias model, proportional to <f>, i.e., S g = bgS + 6^,/nl0 + with all 
the consequences this implies for clustering statistics, e.g., P gg (k) = bgPss (k) + 2bsb t f,fm J P4>s (k) + 
qq ~, ^/nl-FW W + ■•■ • This result is consistent with calculations based on a model for dark matter 



halo clustering, showing that the form is quite general, not requiring assumptions about peaks, or 
the formation or existence of halos. The halo model plays the same role it does in the usual bias 
| picture, giving a prediction for b^ for galaxies known to sit in a certain type of halo. Previous 

projections for future constraints based on this effect have been very conservative - there is enough 
volume at z < 2 to measure /nl to ~ ±1, with much more volume at higher z. As a prelude to 
the bias calculation, I point out that the beyond-linear (in <f>) corrections to the power spectrum 
of mass-density perturbations are naively infinite, so it is dangerous to assume they are negligible; 
qq ■ however, the infinite part can be removed by a renormalization of the fluctuation amplitude, with 

the residual fc-dependent corrections negligible for models allowed by current constraints. 
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I. INTRODUCTION 



Inflation [l], 0, 0> H| has been tested successfully mainly through measurements of the power spectrum of primordial 
density perturbations 0, @] . These perturbations are expected to be nearly Gaussian and nearly scale invariant 
HI, S EH EH ■ Testing the Gaussianity of the perturbations with increasing accuracy will be a major goal of future 
work aimed at distinguishing different models (see [l2| for a review of non-Gaussianity from inflation). The simplest 
models of inflation pre dict that non-Gaussianity will be undetectably small [IH, EH, EH, EH EH Il7l llH . but multifield 
models 119 . 12CI l2ll . 122 . l23l . |2J], models where non-Gaussianity is generated during reheating [25|, [2fJ| or preheating 
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22,123,1301, bouncing/ekpyrotic/cyclic models [31l, I32J, 133J, l3J] , or inflation models based on nonlocal field theory 
36(| can predict levels of non-Gaussianity near the present detection limits (/nl ~ 100, as defined below). While 



there have been some hints of non-Gaussianity in the cosmic microwave backgroun d ( CMB) [371 . |38| |. the general 
consensus seems to be that nothing convincingly primordial has been detected [a, HH, E(| . 

Recently, [4l|, |42|, |43[ showed that there should be a distinctive signature of non-Gaussianity in the large-scale 
power spectrum of dark matter halos, observable as galaxies, in the local model for non-Gaussianity with curvature 
perturbations proportional to 



* = ^ + /nl(0 2 -(0 2 ))+W 3 + ... (1) 

[l2l [44], |4f| [4|| , where I include the 3rd order term which is necessary to compute the power spectrum to 4th order 
in the small, Gaussian, perturbations <f>. This form is only a special case of non-Gaussianity, but it is simple, and 
calculations using it should point the way to generalizations. The idea of [4l|, to use the power spectrum, was a 
departure from the previously standard approach of studying non-Gaussianity in lar ge-s cale structure (LSS) using the 
bispectrum [13, SI, SI, H3, Hll, H3, HI, H3,1HB| or other higher-order statistics [H, 113,111 HI, Hfll, . It was somewhat 
surprising that this signal would be detectable, because the nonlinear corrections to the power spectrum in this model 
were generally assumed to be negligibly small (because fluctuations in (j> are of order 10~ 5 ). Reference [62[ verified and 



further explored this idea and applied it to real galaxy and quasar data sets. They found observational constraints 
-1 (-23) < / NL < +70 (+86), at 95% (99.7%) confidence, by combining LSS and CMB data. The projected rms 
error from the Planck satellite, measuring the CMB anisotropy, is cr/ NL ~ 5 [H, [63[. 
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The calculations of [4j], [42|, [43|, (64] were all based on models for the clustering of dark matter halos. The pur pos e 
of this paper is to investigate the result using a different approach - the renormalized perturbative bias model of [6J] . 
Instead of focusing on dark matter halos, this approach starts by assuming that the tracer density field is a completely 
general, unknown, function of the local mass-density field and Taylor expanding this function in the mass-density 
fluctuations, leading to the form [65j] 

S g (8) = C1 S + \c 2 (S 2 - a 2 ) + ^c 3 S 3 + e + (St) , (2) 

where 8 is the mass-density perturbation, a 2 = (<5 2 ), and 8 g is the tracer density per turbation (I will usually refer 
to the tracer as galaxies, although there are many other possibilities like quasars 62 1, the Lyq forest [66|, [63, HI], 
galaxy cluster/Sunyaev-Zel'dovich effect measurements [6!|, and possibly future 21cm surveys [zO,[zH). The unknown 
coefficients of the Taylor series have become the bias parameters Cj and a random noise variable e has been added to 
allow for stochasticity in the relation (e.g., from shot noise, at least). This type of bias model has been explored in many 
papers [zl, [z!, [zl, (Za, [zl, [z3, [zl, [zl| ■ The new development in [64| was to show that large higher-order corrections 
in this perturbative approach can be eliminated by redefining the bias parameters to absorb them, motivated by the 
way that masses and coupling constants are redefined to absorb divergent loop corrections in quantum field theory 
[SOI ] . I argued that this renormalization approach is an improvement over the previous approach of defining the galaxy 
density to be a function of a smoothed mass-density field (so the higher-order terms in the Taylor series can be kept 
explicitly small [65], [73j]), because the smoothing needed is so extreme that it directly affects the scales of interest 
[HI, HH, and because, in the unrenormalized approach, the higher-order corrections modify the power on arbitrarily 
large, ideally truly linear, scales [zlj]. Reference [83[ found extremely good agreement between the renormalized model 
and the clustering of galaxies in numerical simulations. Generally, the value of perturbation theory (PT) for describing 
LSS will only increase as observational measurements become more precise, because this increases the range of scales 
where corrections to linear theory are important but are still small enough to be treated perturbatively. The value 
of PT has been further enhanced recentl y b y the introduction of several renormalization methods applicable to the 
mass-power spectrum calculation [H iHM H H, US S HI, ]M, 111 ■ 

To complete the relation between primordial non-Gaussianity and final density, note that the initial density field 
is related to through the transfer function from primordial to late-time linear fluctuations, T(fc), and the Poisson 
equation 

8! (k, a) = \p^k 2 T (k) $ k ee M (k, a) 4> k , (3) 

using the definitions of [55(, where I will use the subscript on 8 to indicate the order in the initial density perturbations, 
not the Gaussian field <f> (I will use Sl to indicate the fully linear, including in </>, density perturbation). The growth 
factor is normalized so D (a) = a in the matter dominated era. The transfer function is time independent and 
normalized by T (k — -> 0) — > 1. These definitions make <&k time independent and close to scale invariant. The 
perturbative density field is written as a Taylor series in 81, i.e., 8 = 81 + 62 + 83 + with Si ~ O [S\j, e.g., 

8 2 (k) = / ^3 8, (q) Si (k - q) jf (q, k - q) . (4) 



(2tt)- 



where 



7(2) n 1 ^ 5 1 ki ■ k 2 /fci /c 2 \ 2 /ki ■ k 2 \ 

Js (kl ' k2) = 7 + 2^T U + + 7 {-hkV) (5) 

is given by standard LSS perturbation theory (see [94[ for a review). 

The rest of the paper goes as follows: To put the perturbative bias calculation on a firm foundation, in Sec. [TT1 I 
compute the mass power spectrum in the non-Gaussian model. In Sec. IIII1 I calculate the power spectrum of a biased 
tracer. Finally, in Sec. IIV1 I discuss the results. 



II. MASS-DENSITY POWER SPECTRUM 



In this section I compute the power spectrum of mass-density fluctuations in the non-Gaussian model. I note in 
advance that infrared divergences will appear, i.e., sensitivity to arbitrarily large scales, which suggests that a more 
sophisticated calculation than the usual, essentially Newtonian, LSS perturbation theory might be required. The 
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point of this section is simply to show that the usual LSS PT can in fact be pursued self-consistently, without ignoring 
any divergences, imposing arbitrary cutoffs, or integrating into a regime where the calculations are not valid. All 
of the sensitivity to very large (and very small) scales can be absorbed into the overall normalization of the power 
spectrum, so there is no need to worry about exactly where a physical large-scale cutoff might appear, as long as it is 
much larger than t he scale of the observations. For more physical discussions of infrared divergences in inflation, see 
[§E S [93, H [ll, [Tool, [Toil, [lOi]. Note that there are no divergences in standard LSS PT calculations of the mass 
power spectrum starting from Gaussian primordial pertu rbat i ons, for a typical ACDM power spectrum (divergences 
for pure power law linear power spectra are discussed in [1031 . 1104 . Il05j |). 
The power spectrum to 4th order in 5i is 

Pss (k) = P Sl5l (k) + 2 P SlS2 (k) + 2 P SlSa (k) + P S2 s 2 (k) + ... (6) 

The most interesting term turns out to be the first one, 

P SlSl (k) = M 2 (k,a)P^(k). (7) 

with 

(fc) = P^ (fe) + 2 toL^V + /nl^W (fc) + ... (8) 

The first term here is simply the usual lowest order power spectrum, which I assume is given by P^ (k) = Ak Ua ~^. 
The 2nd term is divergent, in a very simple way: 

2 g^hP^ (k) = 6 <7nl (k) (9) 

where 

(2tt) 
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For n s < 1, er^ — A (e" s 1 — A Ils x ) /2ir 2 (1 — n s ) is infrared divergent, where e is the minimum q and A is the 
maximum q (which we would like to take to zero and infinity, respectively). The solution to this divergence is obvious: 
the divergent term just renormalizes the amplitude of the power spectrum, i.e., A — > A + 6 <7nl cj \^ > - 
The third term in Eq. J5} is slightly more subtle: 

— ^P^(q)JW(|k-q|) [4/£ L ^]P 00 (fc) . (11) 

There is again an infrared divergence for n s < 1, but it does not look like a simple k- independent renormalization of 
the amplitude, because, after we absorb the divergent part, there will be a fc-dependent piece left over, i.e., 

f&Ppp (fc) = P 00 (k) + 2f 2 L [ -^3 [P U (q) P 00 (|k - q|) - P n (q) P^ (k) - P 00 (k) P n (|k - q|)] , 

J (27T) 

(12) 

where I have written the two equivalent negative terms in the integral to emphasize that half of the divergence comes 
from q — > and half from q — > k. 

Subtracting the divergent part has now introduced a new source of trouble, unfortunately, in that the subtracted 
terms very nearly diverge as q — » 00, e.g., for n s = 0.96 a cutoff A > 10 50 fc is required for convergence to 1% accuracy. 
The solution is to pull out another factor renormalizing the amplitude, as follows: 



/^ L P 0202 (fc) = P 00 (k) \Af^a' H + 2f 2 L 

d 3 q 



rf 3 q 
(2^f 



p<m> (q) P<t>4> -q|) p , , p 



(13) 



+ 2/^P^ (k) J 



(2irf 



P00(g)P00(|k-q|) P00(g)P00(|k^ -q|) ' 
p — ttt P^0 (|k - q|) „ — 77-s 1" ^00 ll k * _ ql) 



I have simply added and subtracted P^^, (fc) times some constant piece, with the constant chosen to cancel the 
correction term at some arbitrary scale fc*. Note that we could drop the terms in the integral that are simply integrals 
over —Pcjjcj, (|k — q|) and P^ (|k* — q|), because they exactly cancel. This would leave us with the form we would have 
found by starting with the second form of renormalization (involving fc*) instead of going through the initial step of 
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subtracting the truly infinite part, i.e., the first renormalization was not strictly necessary. I leave the equation in 
the more complicated looking form because the simpler one would involve cancellation between pieces of the integral 
that diverge at different values of q (i.e., fc and fc*), which makes the integral less straightforward to evaluate. 
The power spectrum of <£> is now simply 

P^{k) = P^(k) (14) 

d 3 q 



+ 2/ 2 L P 00 (fc) J 



(27T) 3 



P^(<z)P^(|k-q|) P^{g)P^{\K-g\) , 

p<p<i> (l k - q|) 5 — 77— \ r P<t><t> U k * _ q| 



p^(/0 WVI ^" p^(fc*) 



where P^ now has the renormalized amplitude, essentially disconnected from the original amplitude. Inserting 
P(j>(j> (k) = Ak ns ~ , we find, at least to a very good approximation 

P»« (k) ~ (fc) + / 2 L ^ 2 J? _ i} (fc 2 — 5 - fc"-- 4 *?-- 1 ) = P^ (fc) (l + (? ^ L i} ^ ) . (15) 

Note that this result is well behaved in the scale-invariant limit, (n s — 1) — > 0, when (fc™ 3-1 — fc™ 3-1 ) / (n s — 1) — » 
ln(A;/fc + ) (this would have been a disaster without the second renormalization, introducing fc*, of the UV near- 
divergence, which becomes a true divergence in this limit). 
Now we can estimate the change in apparent slope, 

An s = din (l + ( ^/! L 1) A(fcm '"2 7r T fc "'" 1) ) / dlnfc * ^niA (fc) , (16) 

where A 2 ^ (fc) = Afc ns ~ 1 /27r 2 , and change in the running of the slope 

Aa s = ^^(n s -l)An s . (17) 
dlnfc 

For a realistic model A 2 ^ (fc) ~ 8 x 10~ 10 so the change in slope is quite small, e.g., for an unrealistically large 
/nl = 10 00, An s ~ 0.0033, which is roughly the expected precision of a measurement of the slope using Planck data 
[l06t Il07| . The bottom line of this calculation is that we can safely ignore the non-Gaussian contribution to Ps^, 
for realistic models. Note that, even if the new contribution to the slope or running had been significant, it would not 
necessarily have been distinguishable from a change in the underlying inflation model, i.e., the bare slope or running. 
It would still need to be kept in mind, however, in order to correctly interpret measurements of the power spectrum 
as constraints on the inflation model. 

The rest of the calculation of the mass power spectrum is less interesting, 

2P SlS2 (fc) - 4/ NL f j(2) (q, k - q) P^ L (q) [2P# L (fc) M (|k - q|) + P„ L (|k - q|) M (fc)] (18) 

J (2ir) 

where P^g L (k) = M (fc) P^ (fc). For realistic power spectra, Ps 1 s 2 (k) is comfortably convergent in both limits, and 
produces a subpercent effect for /nl < 10 5 . The terms 4th order in S\ are of course just the usual nonlinear correction 
to the Gaussian linear theory power spectrum. 

To be clear, the conclusion of this section is that the non-Gaussianity of the initial conditions in the local model does 
not significantly affect the mass power spectrum, verifying what has been previously assumed. Note that [881 ] considers 
a different form of renormalized PT with non-Gaussianity, but the results are not directly comparable because their 
specific example is a different form of non-Gaussianity (local in density). 

III. BIAS 

This section contains the primary calculation of this paper: the power spectrum of biased tracers of the density 
field. I start by assuming the bias takes the form of a local Taylor series in density 

5 g = p a /p a - 1 = c s 5 + i c S 2 (S 2 - a 2 ) + i c 5 3 S 3 + e + (Sf) . (19) 
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The first, and only, new nontrivial term related to non-Gaussianity is the linear-quadratic cross-term 

M(k) f -^P , L (q)P^(|k-q|) + 2P^,(fc) / _^p^(g)M(|k-q|) (20 ) 

J (27TJ J (27T) 



C5CS2P55 2 (k) = 2csC52f Nh 

+ /nl— independent parts 



The second term in the brackets is not well behaved in the limit q — > oo, where it becomes [cscpifNLO'ss] -P<Ml (&) 
(this term comes from the part of S (x) S 2 (y) proportional to 8 (x) 5 (y) <f> (y) \7 2 (j> (y)). It nearly diverges, as discussed 
above, in the same sense as the unsmoothed linear density variance nearly diverges, i.e., the cutoff-dependent limit 
goes like A™ s_1 (at least until one hits t he cold dark matter free streaming scale, which actually happens before the 
integral would otherwise converge [X08j | - the fact that that scale enters this discussion should make it clear why 
something needs to be done about this term). It is not immediately obvious how to deal with this term. Ignoring it 
is not a good option, because it can become large relative to the leading order term, at small fc, i.e., 

W„P** (*) _ ^4/ NI ^M- (*) ~ 29 )(??- )(± )( —JL J ) " . (21) 



c 2 s p 5LSL (k) c s J 66 w Vioo A cs J \iooy V0.01/iMpc 

Taking this result at face value is not an option, both because it is a large correction in what is supposed to be 
a perturbative expansion, and because it is ridiculous in LSS perturbation theory to trust the quantity a 2 s which 
requires integration many orders of magnitude into the nonlinear regime for convergence. We cannot remove this term 
by renormalizing the linear bias parameter, because the k dependence of the linear term, proportional to Ps L s L (^)j 
is not consistent with the k dependence of this term, proportional to P$s L (k) (we could always cancel the value of 
this term at one point fc*, but the differing k dependence would mean that the term would become large again very 
quickly as we went away from fc*). 

The solution seems to be to introduce a new term in the bias formula c^,/nl<^! i-e., 

S g = Pg/pg - 1 = Cs 5 + C<p f NL (j) + i c S 2 (S 2 - a 2 ) + ^ c s i 5 3 + e + O (Sf) . (22) 

Then the power spectrum will contain the term 2csc r j,f^P l f ) s L (k), which is needed to absorb the divergence, i.e., the 
bare value of (which might be zero) will be renormalized to b,p — + 2c$2 a 2 s . 

Inevitably, we now also obtain a contribution to the power spectrum c^f^P^ (fc), i.e., the linear theory power 
spectrum is now 

Pgg (*) - c 2 P 5l s l (k) + 2c s c^ L P^ L (*) + c 2 f 2 L P M (k) + ... (23) 

Note that the decision to write c^/nl instead of including the /nl factor in is suggestive but ultimately arbitrary, 
since can in principle depend on /nl- The divergent correction to is of generally the same order of magnitude, 
i.e., proportional to ass, as the correction to eg introduced by [64] . which gives us some reason to expect bs and to 
be of similar order (I use b in place of c to indicate the renormalized parameter). Figure [T] shows an example of this 
power spectrum. 

Equation (|2"3")l is really the bottom line result of the paper, but for completeness I compute all the higher-order 
terms arising from the new term 

f eft 

2c /nl [^P^ (k) + ^fP^ {k)\ = 4c c 52 /2 L P^ (fc) J (?) M (|k - q|) (24) 

Note that each piece here has a near-divergent part, however, each of these is neatly canceled by the renormalizations 
we are already doing. The first term has the same form as the ones we dealt with in Eq. ([21]) . so one can see that 
the redefinition 6^ = + 20,52 a 2 s will cancel the divergent part of the first term in Eq. (I24p . Similarly, the original 



renormalization of cs in [64|, bs = cs + (§72,52 + 2|L) a 2 s , cancels the rest of Eq. ([24)). when bsS is cross multiplied 
with c^/nl^- Things work out so neatly simply because P^x (k) — M^ 1 (fc) Psx (k) (where X could be anything), 
so there is no nontrivial new behavior when calculating cross terms with <p instead of 5. The cross terms between the 
new linear <j> term and higher-order gravitational corrections to 5 are negligible because the latter only contribute at 
small scales, while <p is only significant at large scales (note, however, that if one were trying to exploit the very small 
non-Gaussian signal that is present at relatively high fc, these terms might come into play). 
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k [h/Mpc] 

FIG. 1: Real-space power spectrum for /nl = 5, at z = 0.75, of galaxies with bs = 2, b$ = 3 (black, solid line). Also shown are 
the components of the power spectrum, bgPss(k) (blue, dotted line), 2fjqi J b r / 1 bsP < f,s(k) (green, short-dashed), and f^t,bj,P^(j,(k) 
(red, long-dashed line). The arrow indicates the approximate minimum k probed by an all-sky survey out to z — 1. 

This linear-quadratic (S5 2 ) cross term is the only new interesting one because the linear-linear term is of course just 
the usual linear bias times the nonlinear power spectrum (shown to be unaffected by non-Gaussianity in the previous 
section), while the quadratic-quadratic and linear-cubic terms are already 4th order in <p, without the non-Gaussian 
part. 

After the divergent pieces are subtracted from Eqs. (|2T]) and ([23]). the fc-dependent residuals are tiny for any 
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reasonable /nl, so I drop them. The final power spectrum is thus 




(25) 



where the last three terms are those found in 64], unrelated to primordial non-Gaussianity. 



IV. CONCLUSIONS AND DISCUSSION 



In Sec. HH we found that the corrections to the mass power spectrum arising from the non-Gaussian initial conditions 
are naively infinite, but the divergence can be removed by a renormalization of the amplitude of perturbations, 
after which there is no significant effect of non-Gaussianity for realistic models. It is interesting to note that a 
sufficiently large numerical simulation would have to deal with this issue explicitly, because the most natural method 
of implementing initial conditions effectively implements the unrenormalized version (4lj |; however, the dynamical 
range of typical simulations is small enough that this effect is unlikely to be noticeable (the box size and finite 
resolution provide large and small-scale cutoffs). The issues discussed in Sec. [TT] should be revisited when considering 
other forms of non-Gaussianity. 

In Sec. Mil I computed the power spectrum of biased tracers of mass density and found a large higher-order 
perturbative correction that can only be removed by renormalization if we add a term to the standard linear bias 
model, proportional to cf>, producing the model 



This new term will propagate into the computation of all statistics, not just the power spectrum. Note that <p should 
not be confused with the Newtonian potential at the time of observation - it is defined in terms of S by Eq. ([3]). 

One might ask at this point: "once there is a linear term, why did we not add higher-order terms in (/>?" The 
short answer is simply that we did not need to. We added the minimal term necessary to produce a well-behaved 
calculation. A calculation to even higher order would no doubt require the addition of higher order counter terms, 
and guide their form. It seems likely that the higher-order terms would be completely negligible, like the corrections 
to the mass power spectrum, but, as we have learned, one should really compute them instead of assuming. 

Measuring /nl accurately depends on predicting b^ for an observable in just the same way as measuring the 
amplitude of the matter power spectrum depends on a prediction of the usual bias bs- By comparison to the results 
of ili, we can determine the value of b$ predicted by the halo model. As emphasized by [62| . however, it 
should be kept in mind that this kind of prediction makes some assumptions that are not guaranteed to be true, e.g., 
the galaxy population in a halo might depend on something other than halo mass alone, like merging history, which 
affects the clustering of the population. Halo model predictions will give a good estimate for the size of effect we 
expect for a given /nl and type of galaxy, but if we are fortunate enough to make a detection it will be very difficult 
to measure /nl to high precision using only this power spectrum effect. With that caveat, Eq. (18) of |62fl leads to 



Reference [4J| computed that a very modest future survey, extending only to z — 0.7 (while the bulk of the volume 
of the universe is at higher redshift) can constrain /nl to ~ ±10. It is interesting to know at least roughly how well 
a larger survey can do. For a well-sampled survey (i.e., no shot noise), the signal to noise of a single mode is 



(note that a null test would not include the signal in the denominator, but this makes very little difference, as I will 
show) where 




(26) 



6 = 2 (b s - I) S c = 3.372 {b s - I) . 



(27) 




(28) 




(29) 
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I am assuming T (k) ~ 1 on the relevant scales, as discussed below. I will ignore redshift space, geometric, and 
evolutionary distortions, and any confounding systematic errors or parameter degeneracies. The total signal to noise 
of a survey is approximately 



SV V 



N J At: 2 



dkk 2 



!<*> 



(30) 



where fc m in ~ 2TT/V 1 / 3 (the detection limit is not very sensitive to fc m j n , as long as fc m i n is basically set by the 
survey volume). fc max should be something less than oo, but, for interesting k* the integral converges within the 
linear regime (the integral would not converge if the k dependence of the transfer function was included, however, 
the information gained this way would come from a tiny change in signal on scales where the signal would not be 
uniquely distinguishable from changes in other parameters, and linearity cannot be safely assumed, so a projection 
including it would be unreliable). The bulk of the signal comes from the range 2k+ < k < 10fc*. Consequently, the 
bulk of the signal comes from the /nl</> S cross term, rather than the /n L </> 4> term (cross correlating with an unbiased 
tracer is worse than autocorrelation, but essentially only because of the factor of 2 that is lost in the cross term) . To 
a very good approximation, S/N — 1 when fc* ~ nV^ 1 ^ 3 ~ fc m i n /2 (e.g., this would be ~ 0.8 ttV^ 1 / 3 if I used the 
null test S/N per mode and set k m in based on the diameter of a sphere instead of the side length of a cube, which 
makes it slightly smaller), which gives hhb^/b s D{a) ~ 200(Gpc//i) 2 /^ 2/3 (for Sl m = 0.3). Assuming Eq. (J27J) and 
b s D{a) = 2, we have finally that S/N = 1 when / NL - 125 \(Gpc/hf /V 2 ' 3 } / [2/D(a) - 1]. I plot the value of this 

function for a survey covering all volume up to z in Fig. [2] Note that the assumption bgD(a) = 2 appears to make 
bs unreasonably large at high redshift, however, the value of bs cancels in b^/bs if bg is very large, so the result is not 
actually very sensitive to this assumption. We see that there is plenty of volume in the Universe to measure /nl even 
if it is less than 1, although this will of course take a huge amount of work. Reference [43j found similar results, using 
a variety of slightly different assumptions, which testifies to the robustness of the projection (which ultimately just 
amounts to the statement that the survey should be large enough to resolve the k where th e non -Gaussian-induced 
power becomes roughly equal to the Gaussian power). Twenty-one cm intensity mapping 71, 109] may be a route to 
surveying all of the post-reionization volume. 

Even if one's tracer is weakly biased, so the halo model predicts little signal, it may be possible to make a nonlinear 
transformation of the density field to enhance the bias; however, it is not clear that any such transformation would 
be better than simply measuring the bispectrum. It is interesting to note that the power spectrum measurement 
proposed by [4l[ and discus sed i n this paper can be seen as nature's implementation of the poor-person's bispectrum 
measurement suggested by |llOj |. i.e., the cross-correlation of the field with the square of the field, where here the 
squaring is done for us, as part of the formation of biased structure. The significance of the difference between the 
natural and artificial versions should not be underestimated - exploiting the natural version we can observe the effect 
on large scales even if limited resolution or noise prevents us from probing the small scales that would be required 
to perform a similar transformation artificially, [l lOL llll| found that /nl may be measurable to ~ 0.01 using the 
bispectrum of high redshift 21cm observations and probing down to very small scales. In this paper I have excluded 
the possibility of measuring the signal using small scales (by assuming the transfer function is 1 in the signal-to-noise 
ratio calculation), but it may be useful, especially at high redshift, to explore the possibility of using smaller scales 
and combining the bispectrum and power spectrum. 

I thank Niayesh Afshordi and Latham Boyle for helpful discussions. 
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